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Abstract. We solve the problem of discrete translocation of a polymer through 
a pore, driven by the irreversible, random sequential adsorption of particles on 
one side of the pore. Although the kinetics of the wall motion and the deposition 
are coupled, we find the exact steady-state distribution for the gap between the 
wall and the nearest deposited particle. This result enables us to construct the 
mean translocation velocity demonstrating that translocation is faster when the 
adsorbing particles are smaller. Monte-Carlo simulations also show that smaller 
particles gives less dispersion in the ratcheted motion. We also define and compare 
the relative efficiencies of ratcheting by deposition of particles with different sizes 
and we describe an associated "zone-refinement" process. 



PACS numbers: 05.60.-k,87.16.Ac,05.10.Ln 



1. Introduction 

Polymer translocation through membrane nanopores is a common process in living 
cells. The transport of proteins and nucleic acids in and out of organelles serve a 
variety of control, signaling, and error correction functions [U El [3j . Recent advances 
in polymer manipulation at the nanoscale level have also sparked interest in pore 
translocation as a new tool in genetic sequencing, structure determination, and drug 
delivery 4, 5, 8j. Active polymer transport through pores requires driving forces which 
often provided by "chaperone" proteins that bind to the polymer on one side of the 
membrane. The proteins are larger than the pore, and once bound, create a barrier 
blocking backward polymer fluctuations. This ratcheting process eventually drives 
the entire polymer through the pore. Another known translocation mechanism is by 
"power-stroke" [3J [7J, where chaperones that are deposited close to the membrane 
are subject to conformational changes. These induce a strain that is relieved only by 
direct pulling of the protein through the pore, similar to the driving mechanisms of 
motors such as myosin and kinesin 8 . In post-translational protein translocation, 
both Brownian ratcheting [9l [10] and power stroke [11] pulling, mediated by Hsp-70 
ATPases, have been proposed. Both models exhibit qualitatively similar behavior and 
cannot be distinguished by experimental data [3J. However, translocation by power 
stroke is molecularly more complex, its modeling requires additional parameters [3j [12] 
and its effects arise only at extremely high binding protein densities [7J. Thus, we will 
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Figure 1. Random sequential adsorption of chaperone particles of k lattice sites, 
with deposition rate r. In the polymer reference frame the fluctuating wall has 
intrinsic hopping rates p, q and is ratcheted by the irreversibly bound particles. 
In the case of Hsp70 proteins binding to nucleic acids, k ~ 6. 



only consider Brownian ratcheting as the dominant translocation driving force, and 
neglect power stroke mechanisms. 

The problem of translocating polymers driven by Brownian ratchets has been 
considered by many authors. Analytical progress, however, has been possible only 
for certain limiting cases where restrictions are placed on the binding kinetics. For 
instance, continuum Fokker-Planck models have been solved [H [H GH E] only in the 
limit of perfect ratchets where particles are forced to deposit next to the membrane, or 
in the limit of rapid particle equilibration (Langmuir kinetics). While the continuum 
approach is justifiable in the limit of large chaperone particles that occupy k 1 
lattice sites, in many applications this limit may be unrealistic. For instance, 
chaperones of the Hsp70 family, commonly employed in polymer translocation across 
the endoplasmic reticulum, are approximately 2nm in size. If we assume polymers of 
nucleic acids with interbase distances of ~ 0.36nm, and that the polymer diffuses one 
base pair at a time, typical Hsp70 class chaperones would be described by binding 
particles with k ~ 6. For the translocation of such structures, the discrete approach is 
more pertinent. Random chaperone particle deposition was recently studied using a 
discrete master equation, and including particle detachment and diffusion. However, 
results were derived only in the limit of rapid equilibration, either of the binding 
particles or of the fluctuating polymer [13j . In this paper we report an exact steady- 
state solution of the discrete translocation process in the irreversible deposition limit. 
In our model, the only constraint is that the binding particles do not overlap. No other 
approximations are made. As in previous work we will only consider stiff polymers 
that do not contribute conformational entropy [21 [?1 [13] . The dependence of the mean 
velocity on the size of the deposited particles is explicitly computed. Our result for 
the mean translocation velocity is verified using Monte-Carlo simulations. Simulation 
results for the dispersion of the translocation is also discussed. 

2. Solution of Random Sequential Adsorption Ratchet 

Consider particles of integer length k irreversibly binding to an infinitely long one- 
dimensional lattice representing the translocating polymer. The fluctuating polymer 
is assumed to jump one unit to the right or left with rates q,p respectively. In the 
reference frame of the polymer, as shown in Fig.[l] the membrane wall hops forward 
and backward with rates p, g, respectively. Particle deposition occurs at rate r only if 
there are at least k open sites between the wall and the nearest deposited particle. The 
dynamics of the wall is closely related to that of its nearest gap, since wall fluctuations 
are allowed only if there are enough empty sites for the wall to perform its random 
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motion. Ratcheting occurs when a gap is large enough for a particle to irreversibly 
deposit, preventing the wall from sliding backwards. After deposition, particles cannot 
diffuse nor detach. 

The master equation for the probability density P m ,n(t) for a wall to be at position 
n G (—00, 00) and for the gap closest to it to be of length m € [0, 00) can be derived 
in analogy with problems in Random Sequential Adsorption (RS A) [14"t [T5] . While in 
RSA all gaps are equivalent and one is concerned with the probability of finding a gap 
of length m anywhere along the infinite lattice, here, we are interested in the random 
deposition of particles in the single gap between the wall and its nearest particle, as 
shown in Fig. [I] The time evolution of P m ,n(t) obeys 

00 

P m ,n = pP m -l,n-l + qPrn+l,n+l ~ [P + Q + K m -k + l)H m _ k ] P rn . n + r ^Pj,n, (1) 

j—ra+k 

where the Heaviside function H m -k = 1 for m > fc and zero otherwise. For m = 0, 
Po,n(t) = qPi,n+i— pPo,n+rJ2'jLk Pj,n- The terms with p, q represent the hopping of a 
free, non-interacting wall, which by themselves would lead to a wall drift proportional 
to p— q. The other terms in Eq.[T] describe RSA dynamics. A gap of length m can be 
produced by the deposition of a fc-mer in a gap of arbitrary length j > m+fc. Although 
there are j — k + 1 ways, each with rate r, of depositing a fc-mer in such a gap, only 
one of these choices will lead to the creation of a gap of length m. Similarly, a gap of 
length m can be destroyed by the deposition of any particle of length k within it, a 
process which occurs in (m — fc+ 1) ways for an overall destruction rate of r(m — k + 1) . 
For to < k, no deposition- mediated gap destruction occurs because a fc-mer cannot 
fit into such a small gap. Upon summing P m ,n(t) over all values m € [0, oo) to define 
Q n (t) = ^2 m P n ,m(t), the terms pertaining to the RSA process cancel exactly and 

Qn = P [Qn-1 ~Qn}+q [Q' n+l - qQ' n ] , (2) 

where Q' n = Q n — Po,n is the conditional probability that the wall is at position n and 
that the site preceding it is empty. Upon multiplying Eq. [2] by n and summing over 
the infinite lattice, 

^f±=p-q±Q' j =p- q (*). (3) 

j=-oo 

where (a) = X)^-oo Q'j ^ s the probability that the site immediately preceding the 
wall is empty. This definition implies that (a), the realization-averaged value of the 
random vacancy variable a in the frame of the wall, is also the probability for a gap of 
nonzero length to exist between the wall and the nearest particle. Provided Q'j reaches 
its steady-state distribution in finite time, Eq. [3] defines the steady-state mean wall 
velocity 

v=p-q{a). (4) 

Note that the dependence of the average velocity v on the deposition rate r 
resides completely in the term (a). What remains is to find an explicit steady- 
state expression for (a) = YljQ'j- To this end, we sum Eq.[Tl this time over 
both wall positions 77 G (—00,00) and over gap lengths m' > m obtaining the 
cumulative probability distribution for the first gap to be of length m or larger 
R m = Sm'=m S n= - oo P m,n(t -> oo). We may now recognize that i?i = (a), since 
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the probability for a gap of nonzero length to exist adjacent to the wall is equivalent 
to the probability that gap is of any length except zero. Thus, in order to find exact 
expressions for the velocity in Eq.[5]we need to find R\. Performing the sums over n 
and ml in the steady-state limit of Eq.[TJ we obtain the recursion relation for R m , 

m+fe — 1 

[p + q + r(m -k + l)H m -k+i] Rm = qR m +i + pR m -i - r^^Rj, (5) 

j—m,ax{k,m+l} 

along with the normalization condition Rq = 1. Eq.[S] is solved by introducing the 
z-transform, 

oo 

G(z) = V z m R m , (6) 



-OO 



and its inverse 

Jc zm+ 2ni 

given by the Cauchy integral encircling the origin [17 . Although physically m € [0, oo), 
the z-transform is a sum over all integers, extending the definition to R m <o- in order 
to z-transform Eq.[5]we set H m ^k+i to unity regardless of m. In this case, we find a 
first order differential equation for G(z) which is solved by 

G(z) = -^ TI exp\pz (8) 

For notational simplicity, Eqs. [51 ll3l are expressed with the rates p and q normalized by 
the deposition rate r. The R m arising from inverting Eq.[5]are valid only for m > k — 2 
since the smallest m for which H m -k+i = 1 is m = k — 1 and since, for m = k — 1, 
Eq. [5] contains Rk-2- For k — 1,2 the last term in the exponent of Eq. [8] vanishes and 
the above generating function is valid for all values of m. The Cauchy integral formula 
thus yields R m for all values of m > 0. The integration constant Go can be fixed by 
directly applying the condition Rq = 1. For k = 1 we find 

R fp\ m/2 J P + q+m (2y/pq) , 1 , Q S 
R m = - —j TTT^T - ' k = 1 - 9 

where J v is the Bessel function of the first kind of order v. For dimers, G(z) is identical 
to the k — 1 case except for the substitution q — * q — 1 in Eq.[8] For q > 1 the form 
of R m is the same as in Eq.|Hl with the same substitution in q. If we define Z y = J v 
for q > 1 and Z v = l v for q < 1, where I v is the modified Bessel function of order v, 
for k = 2 we obtain 



wn-x(2ya^ii) fe=2 . (10) 

For g = 1 in the k = 2 case, the direct inverse z-transform gives R m = p m T(p + 
1)/T(p + m + 1), where T is the Gamma function. 

Next, consider the case of larger particles k > 3. The R m >k~2 arising from Eq.[5] 
must now be coupled with explicit solutions for R m <k-3 from Eq.[5l where H rn -k+i = 
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0. Let us derive R m >k-2 from the generating function G(z). If we define as the 
j th term in the Laurent series appearing in Eq.[8j exp ^5Zj=2 z ~ = S^Lo ■> 
we can write i? m for m > k — 2 as 

00 a (,c) W 2 

fl m = G ^ J _ Z a (2y/p\q-l\), fc>3, (11) 

3=0 19 1 

where ct = j+ p + q — k + l + m. These z-transformed solutions are valid only for 
m > k — 2, for which we let R m = GoR m - In order to apply the condition R = 1 
and determine Go, we must connect Eq.[TT]to the m < k — 2 equations in {5J. These 
involve i? m up to m = 2k — 3. Go can thus be determined by 



M 



/ Ri \ 

Rk-3 

GoRk-2 



/ p\ 



(12) 



V G -R 2fc -3 / V J 



where M is the (2k — 3) x (k — 2) transition matrix describing the linear subsystem 
in Eqs.03 



M 



( -{p + q) q o 
p -(p + q) q 

p -(p + q) 







0-1 ... \ 

0-1-1 ... 

-1 -1 -1 ... 

q -1 -1 ... -1 J 



(13) 



Solving Eq.[T2] allows us to determine Go and the exact gap distribution R m for all 
parameters p, q,r,k. 



3. Velocity and Dispersion 

In Figs.[Ua-c) we plot the average wall velocity v = p — qR\ as functions of the 
backward hopping rate q, for various particle lengths k and forward hopping rates p. 
Ratcheting from particle deposition is stronger when q/p is large, the wall motion is 
biased towards the deposited particles, and R\ is small. For fixed kinetic parameters, 
particles of smaller size k yield faster translocation and smaller variance. Smaller 
particles are more effective at translocation due to their enhanced insertion rate 
into the fluctuating gap despite taking smaller steps than larger particles. The 
relative difference of the velocity and dispersion among different particles sizes is most 
pronounced in the strongly ratcheting regime where q/p is large. 

Now consider the dispersion of the ratcheted wall. Upon multiplying Eq. H]by n 2 
and summing over all integers n, we find 



d 

-(n 2 ) = p + q(a) + 2p(n) - 2q ^ nQ' n . 

n=—oo 



(14) 
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Figure 2. (a-c) Exact steady-state velocities as functions of q for r = 1, 
p = 0.1,1,10 at different particle sizes k = 1,2,10. Departure from a simple 
biased random walk occurs for large q/p. The irreversible deposition of particles 
ensures v > 0. For each set of values p, q, the maximal velocity arises for the 
smallest possible particles, when k = 1. 



For this calculation, we cannot simply use (n) = vt as implied by Eq. [4] Although one 
expects the realization-averaged (a) to exponentially decay to its steady-state value, 
an initial transient exists before the distributions reach steady-state. Incorporating 
an "integration constant" arising from this transient, and integrating Eq. [3] (with 
n(0) = 0), we define 



(n(t))= p t-q (a)dt = vt + n . (15) 
Jo 

Similarly, we define 

(n) = ^oc -Tzr = =vt + n . (16) 



The two different integration constants no and n' that embody the initial transients 
do not affect the determination of the steady-state velocity v — d(n(t)) /dt, but 
interestingly, affects the wall dispersion. Upon computing (n 2 ) — (n) 2 = 2Dt, we 
find 

2D=p+q(a)[l-2(n' -no)}. (17) 

The dispersion D retains memory of the transients during which the averaged velocities 
have not yet reached v. This contribution to the dispersion is embodied in the n' — no 
term in Eq. 1171 The term n' — uq is always positive because it takes longer for a wall 
with a gap to reach terminal velocity than an unrestricted wall, leading to a larger 
intercept n' . Therefore p + q(o~) is an upper bound for 2D that is accurate for k = 1 
in the r — > oo limit where (a) — * 0. 

In Figs. [3l we show both the mean velocity v and dispersion D of the wall, 
this time as functions of r, with fixed p = q = 1. As r is increased, not only 
does the mean velocity v increase, but so does the dispersion D. Although the 
wall is "pushed" harder by the rapidly deposited particles, its typical displacement 
also increases to slightly overcompensate the sharpening effect of imposing reflecting 
boundary conditions at each deposited particle. In the r — > limit, we expect the 
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Figure 3. Velocities and dispersions as a function of deposition rate r for fixed 
p = q = 1 and k = 1,2, 3. (a) The mean velocity v for k = 1,2, 3. Both Monte- 
Carlo results and the exact solution are shown together, (b) Monte-Carlo results 
for the dispersion D. For r>j),i) ) the limiting values (from Eq. [55) for D are 
1/2, 14/27, and 13/24 for k = 1, 2, and 3, respectively. 



dispersion to approach that of a diffusing particle with reflecting boundary conditions 
on one side: D = p(l — 2/tt) w 0.3634p for q = p. The Monte-Carlo value for 
p = q = k = 1 and r = 10~ 3 gives D = 0.388 in good agreement with the expected 
result. 

The exact result for the mean velocity given by Eq. [4]and R\ from Eq. [9lor[T0lcan 
be simplified in the slow deposition rate limit r — > corresponding to low chaperone 
concentrations. For example, for k = 1 and q/r = p/r — > oo, asymptotic expressions 
for Ri can be found using the recursion relation 

Ju+l(Z) = -J V {Z) (18) 

z oz 
and the asymptotic limit [16] 

^ J -^^ 3^r(2/3),V3 ' ( 19 ) 

giving 

Ji+2 P /r(2p/r) d\nJ v {z) 



3 1 / 3 r(2/3) fr\ 1/3 , , s 

1 )' ' - as p/r -> oo. (20) 



z—2p/r 



r(i/3) 



•hp/ri^p/r) dz 
Thus, we find 

V =Ml-i?l)~3V3^ r l/3 p2 /3 (21) 

in the r/p—t§ limit. This dimensional result implies that at small chaperone 
concentrations, the mean translocation velocity is proportional to the cube root of the 
chaperone concentration. This nonanalytic limit can also be understood physically by 
considering the typical time t* between succesful chaperone depositions approximated 

by 

1 
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where x(t*) ~ \/pt* is the typical distance between deposited chaperones. Upon 
solving for t* and estimating the velocity as v ~ x(t*)/t* ~ y/p/t*, we find 



r 



1/3^2/3 



(23) 



Our results can also be simplified in the limit of infinitely fast particle deposition, 
where a particle is deposited as soon as the first gap reaches a length of k lattice sites. 
The dynamics then becomes that of the so called burnt-bridge model [THEE [III [21]. In 
the one-dimensional burnt-bridge model, certain links (bridges) separated by k lattice 
sites can be crossed only once by the walk, which then generates a biased motion. The 
speed and the diffusion coefficient can be calculated by either using the continuous 
time version of the discrete time method in [T2] , or by using the general results in [22] , 
as has been done in [21j for an unbiased (p — q) random walk. For r 3> p,q, we obtain 

pk(a- l) 2 



a{a k - I) + k{\ - a) 

(24) 

_pk 2 (a-l) 2 [a 2k+2 +4a k+1 (k + l-ka) + k(l-a 2 )-a(a + 4)] V ' 

2 [a(a k - 1) + fc(l - a)} 3 
where a = q/p. These expressions simplify further in the symmetric case p = q: 

2p 

v = ~, 

k + 1 

(25) 

2 k 2 + k + l 

In the large k limit, appropriate for large binding proteins such as single-stranded 
binding proteins (SSB) with k « 60, the average velocity and dispersion also take on 
simple limiting forms. For a = q/p > 1, the large k limit of v and D given in Eqs. [24] 
are 

pk(a-l) 2 pk 2 (a-l) 2 
V = aW ' °= ■ ( 26 ) 

Since q > p, the drift that tends to close the nearest gap prevents insertion of particles 
even if r is large. Thus, both v and D become exponentially small for large k. If p > q, 
the drift tends to open gaps. However, since very large gaps need to be opened to 
allow insertion of a large fc-site particle, the mean velocity and dispersion approaches 
that of a freely diffusing particle: v = p — q, 2D = p + q. 



4. Annealing and Zone-Refinement 

Finally, we consider the lattice coverage far from the wall, in the long time limit. 
The number of particles adsorbed per length of translocation may be relevant for 
considerations about energetics and macromolecular cost. Monomers will cover the 
translocating polymer behind the wall entirely. Since deposition has been assumed 
irreversible, the deposition of each particle is associated with a large energy cost, 
regardless of size. In this case, monomer deposition may be more costly than 
deposition of larger particles, for which the same energy loss leads to a larger coverage. 
Particles of length k > 1 also allow the presence of empty gaps at saturation, further 
minimizing the number of deposited particles. 
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To quantify a particle cost associated with translocation, consider the infinite- 
time coverage 9 < 1 representing the fraction of filled lattice sites far behind the 
moving wall. For monomers, every site will eventually be filled and 9(k = 1) = 1. 
For k > 1, we compute 9{k) by considering the deposition of a particle into the first 
gap nearest the wall, splitting it in two. One of these daughter gaps becomes the 
new "first" gap closest to the wall, while the other one is now an "interior" gap. 
Particles will continue to deposit into the interior gaps, following rules of deposition 
into a finite length segment [22 [24] , until these gaps reach a fixed coverage. For long 
times we can calculate the total coverage by summing the saturated coverage of all 
the interior gaps. For one particular realization of the random sequential ratchet, the 
creation rate of interior gaps of length m from a first gap of length m' by particle 
deposition obeys N m (t) = rH m r^ m ^k, where N m is the number of interior gaps of 
length m. The Heaviside function prevents the creation of a new interior gap if the 
original first gap is not long enough to accommodate the particle and the new gap. 
Ensemble averaging leads to (N m ) — rR m+ k, where R m +k is the probability for the 
first gap to be larger than m + k and (N m ) is the ensemble average. At long times, 
R m +k reaches steady-state and the generation of interior gaps grows linearly in time 



Interior gaps created by this process are themselves filled by further particle 
deposition. At saturation, an interior gap of initial length m reaches coverage 9 m , a 
quantity that can be calculated by standard RS A techniques [23] [24] . Each interior 
gap is bound one each end by a particle of length k. The total initial length of one 
of these interior gaps and an associated end particle is m + k. After further particle 
deposition into this gap, the number of covered sites reaches m9 rn + k. Upon weighting 
over the ensemble-averaged gap length distribution (N m ), the coverage can thus be 
expressed as 



In Fig. H] we plot 9(k = 3) as a function of q, for p — 0.1, 1, 10. For small q and large 
p, the wall moves forward rapidly, largely independent of deposition. The depositing 
particles rarely interact with the wall and the coverage approaches that of standard 
RSA on an infinite lattice, 9 us A [IS]- If q is large and p is small, the wall stays close 
to the nearest particle, occasionally leaving gaps only slightly larger than k in which 
particles can deposit. Thus, the wall slowly sweeps through the lattice, "zone-refining" 
it by slowing down the deposition process and allowing for more complete filling. For 
intermediate values of p, q, we find that the coverage left by the wall is always between 
the contiguous and RSA limits. 

5. Summary and Conclusions 

In summary, we have found an exact steady-state velocity of the discrete translocation 
problem in the limit of irreversible particle attachment. The exact gap distribution 
R m given by Eqs.[9l[T0l and [TT1 allows us to construct the mean velocity from R\ = (a) 
and Eq.[U Smaller particles yield faster translocation and smaller dispersion, while 
larger particles leave less of the remaining lattice covered. In the protein translocation 
problem, neglecting the dissociation rate as we have done throughout this study is a 
good approximation provided chaperone concentrations >nM. In this limit, the result 
of Elston [3] approaches our velocity given in Eq. [25] For p = q, we find an interesting 



(N m (t)) ~ rtR m+k . 
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T,Z=o( m + k ) R m+k 



(27) 
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Figure 4. Particle coverage 6 behind the wall as a function of q for k = 3, 
r = 1, and various p. For walls that move forward rapidly, p — » oo, the coverage 
corresponds to that of irreversible RSA on an infinite lattice without a wall and 
8 RSA (k = 3) = 3-v/7re _4 (Erfi(2) - Erfi(l))/2 = 0.823653 (dotted line) For p -> 0, 
the wall slowly sweeps across the lattice, allowing more contiguous deposition 
behind it. In this case 0(k,p — ► 0) ~ 1. 

r 1 / 3 dependence of the mean velocity for small r. We also find that the coupling 
of particle deposition to wall dynamics allows the wall to "zone-refine" the lattice, 
producing long time particle coverages Orsa < < 1. 

Our results can be used to guide experimental systems that probe the 
mechanisms of chaperone-assisted translocation. For example, comparison of the 
mean translocation speed v with our exact solution and the dispersion D with 
our Monte-Carlo results, lead to independent values for p and qR\ defined in this 
paper. By measuring v and D for different driving force F oc lng (by tuning e.g., 
a transmembrane potential) one can numerically determine both the particle size k 
and the effective adsorption rate r. If particle detachments occurring at rate rd are 
also considered, a sufficient condition for our solution to the mean velocity v to be 
accurate is rd <C vOasA/k. For larger detachment rates, we expect a stall force, where 
the mean velocity vanishes, and may be negative for large enough F. In this case, 
the ratcheting occurs with pawls that detach, allowing occasional backsliding to the 
second particle. 
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discussions. MD and TC were supported by the NSF through grant DMS-0349195 
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